Conceptual Design of a Micron-Scale Atomic Clock 
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A theoretical proposal for reducing an entire atomic clock to micron dimensions. A phosphorus 
or nitrogen atom is introduced into a fullerene cage. This endohedral fullerene is then coated with 
an insulating shell and a number of them are deposited as a thin layer on a silicon chip. Next to 
this layer a GMR sensor is fabricated which is close to the endohedral fullerenes. This GMR sensor 
measures oscillating magnetic fields on the order of micro-gauss from the nuclear spins varying at the 
frequency of the hyperfine transition (413 MHz frequency). Given the micron scale and simplicity 
of this system only a few transistors are needed to control the waveforms and to perform digital 
clocking. This new form of atomic clock exhibits extremely low power (nano watts), high vibration 
and shock resistance, stability on the order of 10 -9 , and is compatible with MEMS fabrication and 
chip integration. As GMR sensors continue to improve in sensitivity the stability of this form of 
atomic clock will increase proportionately. 
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I. BACKGROUND 



Atomic clocks have been in existence since 1949. The 
basic time keeping element is the hyperfine interaction 
between outer electron(s) and the nuclear spin. Thus 
embodiments typically use atomic hydrogen, the alkali 
metals, and ions that have a single outer orbital elec- 
tron remaining; however hyperfine splitting is present in 
more general atomic and molecular species and can be 
used for time keeping. The hyperfine interaction is ro- 
bust against perturbations from vibration or temperature 
since it only involves the density and spin of the electron 
wave function at the nucleus. High quality atomic clocks 
have time precision of better than 1 part in 10 15 . Atomic 
clocks have been reduced to a few cubic millimeters in 
size. There is a great deal of interest by DARPA and 
the electronics industry in reducing the size, complexity, 
and power requirements of atomic clocks. A precise and 
stable time base that can fit within the packaging and 
power envelope of modern devices will greatly increase 
the efficiency and robustness of mobile computing and 
sensor elements. Spectral bandwidth allocation is funda- 
mentally limited by the accuracy and stability of the ref- 
erence frequency used to define the clock period. One of 
the fundamental limits to reducing the size of electronic 
wireless sensors (motes) is the fact that crystal oscillators 
cannot be shrunk beyond the current package sizeQ. A 
MEMS-scale atomic clock could replace complex clock 
synchronization methods within computers and improve 
arbitration protocols between chips by supplying precise 
local clocks. 
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A. Atomic Clock Principles 

Currently atomic clocks operate by using vapors of Cs 
or Rb coming off an oven at ~80 °C. This results in 
about 1 atom per cubic micron. In some optical clocks 
buffer gas is used to slow the collision rate of the atoms 
with the wall where spin flipping occurs - broadening 
the hyperfine line width. Each atom produces a mag- 
netic field due to the spin of the nucleus - 10 -11 gauss at 
1 micron distance. When the atomic system has been ex- 
cited, the magnetic field oscillates at the hyperfine tran- 
sition frequency: 100 MHz to 20 GHz - depending upon 
the atomic species in use. The best current systems use 
stabilized lasers to interact with the atomic system by 
means of optical transitions. Thus they do not directly 
sense the nuclear magnetic field. This method is costly 
in size and power. 



II. ENDOHEDRAL FULLERENES 

Buckyballs (Ceo and C70, each about 1 nanometer 
in diameter) have been developed that contain a single 
atomic passenger and even simple molecules [7l [l2l Il5l. Ho. 
l2lj |. Figure [T] shows a simple model. For most of these 
endohedral fullerenes a charge transfer of the enclosed 
atom to the fullerene cage occurs resulting in a chemi- 
cal bond and distorted structure. In the case of group 
V atoms the trapped atom is confined by a harmonic- 
like potential to the center of the cage. The atom is 
trapped within the covalent bonds of the Buckyball with 
its outer electrons symmetrically repelled away from the 
walls of the cage. In almost all respects the atom be- 
haves as a 'free' (unbonded) atom, though spatially re- 
stricted to be within the fullerene cage. The distributed 
7T bonding electrons also act as an almost perfect Fara- 
day cage. The inner atom is fully capable of entering 
into magnetic interactions since the Buckyball cage is 
spin neutral (it is desirable to use C 12 to remove nuclear 
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FIG. 1: Model of a Buckyball with an enclosed atom 

spins within the cage walls). Two useful compounds are 
N@C6o and P@Ceo, i.e., nitrogen and phosphorus inside 
Ceo ■ Group V atoms are paramagnetic due to their half- 
filled p-orbitals. Electron Spin Resonance experiments 
have demonstrated that the trapped electrons for group 
V atoms have very long spin relaxation times [18(. 

Group V atoms can be implanted into fullerenes by 
simultaneous ion bombardment and fullerene evapora- 
tion onto a target [l8j. The resulting endohedrals are 
themally[l3| and chemically stable @ at ambient con- 
ditions but do show disintegration between 400K and 
6OOK0. In low-energy implanted samples only a small 
fraction are actually filled, typically 1 in 10,000. It is pos- 
sible to enrich and purify the fraction of filled molecules 
by high-pressure chromatography -8j. Recent work has 
shown that nuclear implantation techniques (using 18 
MeV ions) can convert almost 6% of a thick fullerene 
layer (^40 fxm) into endohedral compounds @. A C 60 
Buckyball with an atom trapped within has a nominal 
diameter of 1 nanometer. 



A. Endohedral Fullerenes Hyperfine Interactions 

Phosphorus comes in one stable isotope: 31 P(100%), 
spin/parity = 1/2+. Its electronic ground state, 4 53, 
is split by the hyperfine interaction. Nitrogen comes 
in two stable isotopes: 14 P(99.6%), spin/parity = 1+ 
and 15 N(0.37%), spin/parity = 1/2-. Both isotopes have 
their ground state split by the hyperfine interaction. For 
the sake of concreteness we focus on phosphorous. Nitro- 
gen may require isotopic separation for use in an atomic 
clock, or, we could accept a double resonance mode of 
operation. 

The contact part of the hyperfine interaction is 

U = 7^b|V>(0)| 2 / • S (1) 
where J is the nuclear spin in units of h ( S is the electon 



spin). Typical values of the hyperfine frequency are in the 
microwave range. 

One important effect is the increase in electron density 
at the nucleus due to confinement inside the fullerene. 
Recent experimental research shows an enhanced electron 
capture decay rate in 'Be when it is encapsulated inside 
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FIG. 2: EPR spectra, [TJ 

Figure [5] shows EPR experimental data[l4| on the 
P@C6o ground state splitting due to the hyperfine inter- 
action. The doublet splitting (A=49.2 G) in the lower 
panel is due to the hyperfine interaction of the electron 
spin with the nuclear spin 1=1/2 of 31 P. The EPR spin- 
magnetic field interaction energy is 

U = ±3[* B B a (2) 

This implies a hyperfine frequency of 413 MHz. This 
is a fairly low hyperfine frequency but reflects that fact 
that the outer electrons have finite angular momentum, 
consequently they have zero spatial overlap at the nu- 
cleus. The hyperfine interaction in paramagnetic atoms 
or ions is entirely determined by the admixture of elec- 
tronic s-orbitals from excited configurations which result 
from atomic interactions such as electron repulsion and 
correlation [l|. 

B. Atomic Transitions 

Given that P@C6o is a paramagnetic atom the most 
straightforward method for measuring the hyperfine en- 
ergy is to polarize the atom in a large static magnetic 
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Hyperfme splitting 



Zeeman splitting 
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FIG. 3: Zeeman/hyperfine energy levels, label = \Seiectron, 
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TABLE I: Nuclear-Electron Perturbation Effects 



Distance [nm] 


Frequency Shift 


1 


10 KHz 


10 


10 Hz 


100 


10 mHz 


1000 


10 juHz 



field. Figure [3] shows the atomic energy levels when the 
electron spin - magnetic field interaction energy exceeds 
the hyperfme energy splitting. In the absence of EPR 
pumping the vast majority of atomic systems will be in 
the lowest energy state. The polarized outer electron 
spins will define a spatial direction for the nuclear spin 
and will define its energy transitions. Usually hyperfme 
transitions inside paramagnetic atoms and ions have a 
poorly defined line width - due to the vast number of 
perturbations acting on the valence electrons, which are 
directly relayed to the nucleus through the hyperfine in- 
teraction. In the case of P@C6o the exceptionally long 
electron spin relaxation times and the nearly ideal con- 
ditions of isolation give us relief from these problems. 



C. Interactions Between Endohedral Fullerenes 

One issue in the use of a packed array of P@C6o 
molecules for an atomic frequency standard is the pertur- 
bation nearby spins apply to each other. Nearby electron 
spins will contribute varying magnetic energy terms to 
the nuclear spin. As each nuclear spin will see a slightly 
different electron spin neighborhood this will broaden the 
hyperfine frequency of the ensemble. The magnitude of 
the shift for each electron-nucleus pair is 
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±3 
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(3) 



Table U shows this interaction for an electron-nuclear 
pair as a function of separation. 

This implies that a close-packed solid of endohedral 
fullerenes (diameter = lnm) will exhibit a fractional fre- 
quency shift on the order of 10 N nearest neighbors- In itself 
this is not a problem as this is a constant shift. The prob- 
lem stems from the variability of the packing. Random- 
ness in the distances to nearby endohedral fullerenes will 
modulate the perturbation. This random perturbation 
will broaden the absorption line of the system resulting 



TABLE II: Density Effect on the Standard Deviation of the 
Frequency Shift(i? cuto ff = lnm) 



Density of fullerenes [m 3 ] 


^/nuclear shift) [ Hz l 


10 17 


5 


10 18 


17 


10 19 


55 


10 20 


172 


10 21 


545 


10 22 


1,721 


10 23 


5,384 


10 24 


16,086 


10 25 


39,288 


10 26 


62,452 



in poor short-time frequency control. Assuming stable 
packing the long term frequency control of the system 
will still be exemplary. 

Since the perturbation term varies as ^a simple solu- 
tion could be to dilute the endohedral fullerenes inside a 
neutral carrier material. The issue now is to estimate the 
effect of random spacing between endohedral fullerenes 
and the corresponding perturbation effects so that the 
net line width of the system meets the short-time fre- 
quency requirement. The main problem is the variance 
in the frequency. Appendix [Al calculates this effect with 
the results in Table ILTl 

The conclusion from this calculation is that the di- 
vergent -^j interaction amplifies the effect of statistical 
fluctuations inside a homogenous mixture. While the ef- 
fects of randomness in fullcrcnc locations limit the ul- 
timate precision of an atomic clock, for less demanding 
uses, e.g., crystal clock replacement, a simple layer of Ceo 
doped with P by nuclear implantation techniques should 
prove an attractive option. 

A random mixture of endohedral fullerenes is not the 
best spatial organization for a frequency standard. The 
best arrangement is a uniform grid with a few 10's of nm 
spacing between occupancy sites. 



D. Classification of Endohedral Fullerenes 

It is possible to separate each endohedral fullcrcnc from 
its neighbors by coating it with a glass shell. Silica 
gel, an inorganic polymer, has a three-dimensional net- 
work and can easily be synthesized via the sol-gel route. 
Fullerenes cannot be incorporated into sol-gel glasses ho- 
mogeneously due to low solubility. This problem can be 
overcome by functionalization of the fullerenes with such 
groups as will form some kind of bond (hydrogen, van der 
Waals, or covalent) with the growing silica net work (l9j. 
The synthesis of silica-fullerene hybrid materials is done 
in the following manner. Silica precursor is mixed with 
alcohol and water. Acid is added as a catalyst. Func- 
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tionalized fullerene is added to this mixture directly or 
by dissolving in toluene. The mixture is allowed to gel 
over a period of time at elevated temperature. Coatings 
of single fullerene molecules can be accomplished by per- 
forming the above steps in a micro-fluidic mixing cham- 
ber where the timed release of materials results in small 
globules of glass-encased single fullerene. 

Fullerene functionalization consists of covalently bond- 
ing side groups mainly consisting of nitrogen or Si-O-R 
(R= Me, Et) or even hydroxyl groups. The molecular 
carbon allotrope readily adds nucleophiles and carbenes 
and participates as the electron-deficient dienophile com- 
ponent in many thermal cycloaddition reactions such as 
the Diels-Alder addition. In all mono-adducts formed 
by 1,2-addition the fullerene preserves the favorable 7r- 
electron system of C6o[5(, however, we can expect a loss 
of the full symmetry in Buckyball shielding - the net ef- 
fect on the hyperfine energy of the system needs to be 
determined by experiment. 
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FIG. 4: Functionalization of a fullerene for sol-gel coating, 

El 




FIG. 5: Functionalized fullerene inside a sol- gel network, [T^] 

So we see that a good arrangement for the endohedral 
fullerenes is to coat them with ~ 20nm diameter glass 
coating. This separates the electron-nuclear spins to the 
point where a narrow line width is possible. It should be 
noted that wrapping polymers around fullerenes is also a 
separation option. 

Another embodiment which avoids the possible prob- 
lems of functionalization is to create simple fullerenes and 



then implant them at low energy into an insulating ma- 
trix. We can also use patterning (e.g., 2D interference 
patterns from hard or soft UV light) to create openings 
in a photoresist layer. Next we implant a low flux of 
fullerenes so as to only deposit at most one fullerene per 
opening. Follow by a thin layer of an insulating material 
and so on. Finally the photoresist layer with the extra 
fullerenes is burned off by an oxygen plasma. There are 
many techniques used in electronics that can create well- 
separated fullerene matrices. 

E. Magnetic Parameters for Encapsulated 
Endohedral Fullerenes 

If we assume a uniform sample of ~ 20nm diameter 
glass spheres with a single fullerene inside we predict an 
internal magnetic field from the nuclear spins of 

-Bnucicus = ^density = 10~ 10 [tesla] = 1/i gauss. (4) 

47T 

At the same time the internal magnetic field from the 
electron spins is 

-^electronic - M-B -^density = 3 x 1CT 7 [tesla] = 3 milligauss. 

(5) 

Because the source of the polarizing field for the nu- 
clear spin is the electron spin coupled through the hy- 
perfine interaction we need to re-express this in terms of 
the effective magnetic field being applied to the nucleus. 
Translating the 413 Mhz hyperfine frequency into a nu- 
clear Zeeman splitting gives an equivalent magnetic field 
of 27 tesla. For the coupling term we have 

^electron — 3/2, 

t/Hyperfine = 4.13xl0 8 [Hz] hi -S= 2.75xlCT 25 I ■ S[Joules], 

„ ^Hyperfine,-, 54[Tesla]S 

-Deflective — * — J, — 9*. (Dj 

UN 6 

It should be noted that the interaction energy repre- 
sented by this pseudo magnetic field is very precisely de- 
fined - it is the hyperfine energy. This expression for 
the interaction will allow us to directly employ the Bloch 
equations below to predict electron spin precession per- 
turbations on nuclear polarization. 

One of the design challenges for the atomic clock is 
to "pin" the electron spins to minimize their precession 
which causes interference with the nuclear signal. Also, 
electron spins are very susceptible to environmental noise 
so we need to ensure that most of the excitation energy 
for the hyperfine defined nuclear spin flips come from 
an external AC drive field directly - not through the 
precessing of the electron spins. 

1. The Bloch Equations 

In our approach there are three spin mechanisms at 
work. The first mechanism is NMR where the DC polar- 
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ization field comes from the spin polarized electrons via 
the hyperfine interaction and nuclear spin flips are driven 
by magnetic resonance with an external AC magnetic 
field. By design this mechanism is on resonance. Bloch's 
equations describe this situation accurately as there are a 
large number of noninteracting spins undergoing preces- 
sion. The second mechanism is induced precession in the 
outer electron spin with an external DC magnetic field 
and the same AC magnetic field as for the NMR mecha- 
nism. With a large Zeeman splitting this precession is far 
off resonance and the electron net spin will only precess, 
never flip. We need to reduce this signal to a small value 
that won't interfere with the nuclear hyperfine signal. 
The third mechanism is nuclear spin precession driven 
by the precessing outer electron spin through the hyper- 
fine interaction. Bloch's equation, using the equivalent 
magnetic field of the hyperfine interaction, works here as 
well. We need to ensure that this mechanism is small 
compared to the conventional magentic excitation. In all 
three cases we can use the Bloch approach. 

Bloch's phenomenlogical equation is an adequate de- 
scription for both the electron and nuclear spin 



— = 7 MAH yJ - 

at T 2 Ti 



M , 
-k , 



(7) 



where k' are the unit vectors of the laboratory 

frame of reference, T\ and T 2 are phcnomenological relax- 
ation times, and the various M's are the magnetization 
vector components. 

In the prototypical magnetic resonance experiment 
a large DC magnetic field, H , is applied along the 
k' direction. An rf magnetic field, Hi, is applied along 
either the i' or j'direction. 

After a series of calculations AbragamQ derives the 
following time-dependent magnetization in the labora- 
tory reference frame 

My = jHxMqJtz (Aw) cos cut, 



/r a (Aw) 



T 2 



7T 1 + Aw 2 T 2 ' 



Alu = LO — U) 

7 = gyromagnetic ratio, 
Mq = the magnetic moment of the sample, 
ujo = J H , 
lo = rf oscillation frequency, 
Hi = magnitude of rf field 



(8) 



This assumes negligible saturation, that is 
7 2 iJ 2 TiT 2 <C 1. For the three spin mechanisms in- 
volved we replicate the Bloch approach, using the 
nuclear parameters in the first case, then the electronic 
parameters, and finally the hyperfine driven precession 



with Hi now coming from the electron spin through the 
hyperfine effective magnetic field. 

Twamley [22[ reports that the electronic relaxation 
times for Group V endohedral fullerenes are Ti ~ Is at 
T ~ 7°K, while T 2 ~ 20/is for all temperatures. He also 
states that nuclear relaxation times should be several or- 
ders of magnitude longer than the electronic relaxation 
times. We will assume 1 second applies for the nuclear 
relaxation parameters at room temperature. 

For the pure nuclear NMR spin mechanism this implies 
that Hi -C 10~ 8 tesla, or 100 micro-gauss. Equation [5] 
implies a resonance when the applied rf frequency equals 
the hyperfine transition frequency. The line width at half 
maximum is Given the other parameters we predict 
an NMR magnetic signal on resonance of about 1% of 
the applied rf magnetic field amplitude. 

For the electron precession effect we use the above 
equation in the regime where the precession frequency is 
far above the rf frequency (for the precession frequency 
to be much larger than 413 MHz the DC magnetic field 
must be much larger than 49 gauss). 

My = 7fliM /r 2 (Auj) cosuit, 



Alo = -ui = - 7 i?o, 
Ti 1 „ 1 



7T 1 + A^ 2 T 2 7T 7 2 /f 2 T 2 



(9) 



From © we predict an electron spin precession mag- 
netic signal of 



3^o 

47T 



1 



^electron = ^ ;p N ~ COSLjt , 



Hi^-fl B N dcnsi ty ^pr— C0S Ujt ' 

4tt ^lH^T 2 



7 



~ h ' 

Hence 



^electron ^ g _M0 hN 1 CQSUjt M()) 

1 uLnsny jj2rp electron ' v ' 



We note that this signal drops as one over the DC 
magnetic field squared. For a 100 gauss H field and 
Hi ~ 10~ 8 tesla we predict an electron precession mag- 
netic field of ~ 10 -18 tesla. This is clearly not a problem 
against the 10~ 10 tesla NMR output signal. 

For our third spin mechanism (the hyperfine pseudo 
magnetic field) we see that an electron spin precession 
signal of ~ 10 -18 tesla versus the ~ 10~ 7 tesla total 
electron magnetic field in the fullerene sample implies 
the fractional electron spin moment precessing about the 
DC magnetic moment is around 3 x 10~ 12 . From Equa- 
tion [5] we predict that perturbation at the nucleus is 
10~ n tesla. This is much smaller than the NMR rf field 
strength. 

Thus to establish a strong and stable electron spin sys- 
tem within which to monitor nuclear spin transitions we 
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need an external DC magnetic field of at least a few hun- 
dred gauss. Pure iron has a magnetization of 1700 gauss 
at room temperature. Thus a simple ferromagnet can 
provide the polarization field. 



III. MAGNETIC SOURCES 

While a coil of wire is a convenient source of uniform 
magnetic fields when dealing with macroscopic dimen- 
sions, this is not true on the micron scale. Fabrication of 
coils in IC's demands the use of vias and complex mate- 
rials fabrication. A better approach is to use Maxwell's 
equations and derive the magnetic field from the displace- 
ment current inside a capacitor. 



A. Dimensional Analysis 

To show that the use of capacitors and typical CMOS 
operating voltages can supply us with an adequate rf 
magnetic field we first perform a crude calculation of the 
rough scale of the phenomena. We assume a micron scale 
for the critical dimensions, 



e ~ 10" n F/m,// ~ 10~ 6 ii"/m,d ■ 



l/zm=10" 6 m 



This implies some characteristic sizes for capacitor and 
inductors 



C ~ 10 _1 'F,L ~ 10~^H 

We estimate the drive current a micron-size inductor 
needs to create a magnetic field sufficient to saturate the 
spin ensemble, 10 _8 tesla. 



$ = BA = LI = 10~ 12 HI 
A ~ 10" 



12 2 

m , 



B=10" 8 tesla. 



Hence I ~ 10~ 8 amps. We now want to determined the 
drive voltage a typical micron scale capacitor must have 
to create this amount of current through the displace- 
ment current given by Maxwell's equations. 



Q = CV, 



I=10" H amps = Q = CV = CuV, 

uj r~ 10 9 rad/ sec 

Hence V ~ lvolt. Thus applying an rf voltage of 
around 1 volt to a typical micron scale capacitor will 
generate the needed 10 _8 tesla magnetic field needed to 
spin polarize the endohedral fullcrcncc atomic standard. 

Given that this calculation is simply dimensional anal- 
ysis we conclude that choices in geometrical ratios and 
materials constants, i.e., permittivity and permeability, 
can change these magnitudes by several orders of magni- 
tude - in both directions. 



1. Example Magnetic Driver 

A concrete example of the above magnetic generator 
is given below. 
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FIG. 6: Balanced capacitors magnetic source 

Here two parallel plate capacitors (the long dimension 
of the plates into the plane of the paper) generate an AC 
magnetic field between them. The source of the magnetic 
field is the displacement current inside the dielectric lay- 
ers, i.e., the time varying electric field. The scale of the 
external field's spatial extent is about the thickness of 
the dielectric spacers. Using two capacitors driven out 
of phase generates a low gradient magnetic field between 
them. This is a two-dimensional version of Hclmholtz 
coils. 



IV. GMR SENSORS 

The GMR effect takes place in heterogeneous magnetic 
systems with two or more ferromagnetic components and 
at least one nonmagnetic component[10]. An example 
is the trilayer permalloy/copper/permalloy system. The 
GMR coefficient for a multilayer system is defined as the 
fractional resistance change between parallel and antipar- 
allel alignment of the adjacent layers. This coefficient can 
be as high as 10% for trilayer systems and more than 20% 
for multilayer systems. 



Zero-field high resistance High-field low resistance 



FIG 7: The Giant Magnetoresistance effect is due to the large 
difference in electrical resistance between two magnetic states 
of a metallic multilayer film 

The occurrence of the GMR effect depends on the abil- 
ity of the applied magnetic to switch the relative orien- 
tation of the magnetic moments back and forth between 
the parallel and antiparallel states. In some multilay- 
ers a quantum-mechanical interlayer exchange coupling 
across Cu or another paramagnetic metal causes a zero- 
field antiparallel alignment which can be overcome by a 
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high applied field. The magnitude of the GMR effect 
can be surprisingly large, up to 80%. However, the fields 
needed to saturate Co/Cu multilayers are too large for 
sensor applications. Other multilayers are designed to 
have an antiparallel state in a limited applied field range 
by alternating ferromagnetic layers (Co and Fe layers in- 
stead of two NiFe layers) with different intrinsic switching 
fields. Thus the behavior of a multi-layer GMR stack can 
be tuned to varying DC magnetic bias conditions. 

Another interesting development is the ballistic mag- 
netoresistance effect, BMR, in ferromagnetic nanocon- 
tacts. The BMR effect arises from nonadiabatic spin 
scattering across very narrow, ~ atomic scale, magnetic 
domain walls trapped at nano-sized constrictions. In 
one study, the observation of a remarkably large room- 
temperature BMR effect in Ni nanocontacts was re- 
ported. The observed BMR values are in excess of 3000% 
and are achieved at low switching fields, less than a few 
hundred oersteds. Q 
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easily rejected if it is directed perpendicular to the multi- 
layer planes of the GMR sensor stack - the various mag- 
netic layers only have substantial magnetic polarization 
in their plane. The driving AC magnetic field is more 
difficult to reject for multi-layer GMR stacks as it will be 
aligned along the planes of the stack, just as is the atomic 
clock magnetic signal. If we choose to detect the NMR 
signal that is perpendicular to the driving field we have 
a chance to selective shield the GMR sensor from the 
AC driving field direction by the use of a slitted Faraday 
shield. 



f — GMR stack 



(Side view) ^\ 

/ 



Copper shield 



Applied field (Oe) 



(Top view) 



FIG. 9: Faraday shielding for GMR directional selectivity 

The use of a slitted copper enclosure greatly reduces 
the GMR sensitivity to magnetic field aligned perpendic- 
ular to the direction of the slits while perserving sensitiv- 
ity to parallel magnetic fields. This holds only for high 
frequency magnetic fields and is due to the induced eddy 
currents that counter the external applied field. 



FIG. 8: Successive BMR loops from a Ni nanocontact showing 
3150% BMR, H 

TMR (tunneling magnetoresistance) and CMR (colos- 
sal magnetoresistance) sensors are also devices that can 
sense magnetic state with increased sensitivity and lower 
power. 

GMR read heads for hard disks operate at gigahertz 
frequencies. This operation frequency is near the ferro- 
magnetic resonance of these devices. One study of these 
devices measured a resonance frequency of 3.60 GHz[2fJ. 
A recent experimental report shows a GMR sensor sen- 
sitivity at room temperature of 350 pT/Hz^fTll]. 



V. APPARATUS 

We now can put together a conceptual design for the 
atomic clock. 

. . GMR sensor i 1 



Drive capacitors 



Fullerenes 



(Side view) 



FIG. 10: Conceptual atomic clock design 



A. GMR Directional Sensitivity 



For atomic clock application it is desirable that the 
GMR sensor reject the magnetic field coming from the 
both the strong DC polarizing magnetic and the AC driv- 
ing field that induces NMR resonance. The DC field is 



A. Atomic Clock Design 

Putting together all the ingredients discussed above 
we have a micron-scale scheme that takes a drive voltage 
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near the hyperfine resonance (driving the capacitors) and 
produces a small (~ 1% of the drive field) signal that is 
detected and amplifed by the GMR sensor into a reason- 
able voltage. Exact device values are highly dependent 
upon materials and fabrication and won't be calculated 
here. In terms of basic physics we expect a fullerene sig- 
nal on resonance (413 MHz with a line width ^ 1 Hertz) 
on the order of 10 _10 tesla and GMR sensors at room 
temperature have achieved a sensitivity of 350 pT/Hz 1 ^ 2 . 
Thus in a 1 Hertz filtered loop we should detect an atomic 
clock resonance with S/N > 1. 



Electronics 



would give us a very well determined perturbation situ- 
ation that can be exploited for accuracy. In the limit 
of true nanotcchnology the ultimate clock is a single 
fullerene with considerable shielding. This should be 
competitive with very good atomic clocks of vastly more 
volume. 



APPENDIX A: VARIANCE OF FULLERENE 
ELECTRON-NUCLEAR INTERACTION 

The idealized variation in frequency shift due to nearby 
endohedral fullerencs 

n = density of fullerenes, 




AGC 



Atomic resonator 



(N) = expected number of neighbors within R{Al) 



nAirr dr, 



(N) = 



(A2) 



(A3) 



Therefore 



FIG. 11: Simple clock circuit 

The above figure shows a simple oscillator design that 
will provide an ouput sine wave locked to the line fre- 
quency of the fullerene clock standard. The AGC circuit 
is designed to provide stable oscillations just below sat- 
uration. Starting from Power On, the intrinsic noise in 
the circuit will provide a nanovolt or so in the 1 Hz reso- 
nance bandwidth. If the amplifier ultimately provides a 
drive signal of around 1 volt with a mild amount of net 
circuit gain, we should have stable full-power operation 
in ~ 100 system time constants, ~ 100 x 1 sec (1/band- 
width). Faster turn on can be provided by strobing the 
input of the resonator with a pulse providing substantial 
fouricr components near the resonance. 



C. Clock Scalability 

The simple scheme we've discussed gives us a micron- 
scale atomic clock with 10 _9 accuracy and a power dissi- 
pation of a nanowatt (10 nW capacitive drive but we can 
use resonant circuits to store the energy) . This will likely 
be adequate for many mobile/sensor net applications but 
not adequate for more demanding situations. What can 
be done? 

First, as GMR sensors improve (BMR, etc.), we can use 
more diluted fullerene stacks to gain a sharper line by a 
cubic factor in separation as we lose an equal amout of 
magnetic signal. A nanoscale-precise placing of fullerenes 



R 



3 
47m 



i(.v)-i = \; — • ( A4 ) 
The expectation value for the nuclear frequency shift 



(/nuclear shift) - 3 M ° M ^ M ™ (K) ■ < ^ ' 



\\-Kr 2 dr = 47m 



t<N) = l ^ 



-dr, (A6) 



ficutoff 



: 47m [log(^Ar)=i) - log(-Rcutoff)] = * n 



(A7) 



<7,i,= variance 



= 47m / (— - k mcan ) r dr 

-'-Rcutoff r 



(A8) 



2 



47m 1 2 3 1 

= ~ 5~ V Trt I" ^mean-R(/V) = l + ~H3 V^) 

It (N) = l -"cutoff 

"^mcan^cutoff _ 6fc moan log[i?(jv) = l] 
+6fc 

mean log[i?cutoff]). (A10) 
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